High-resolution spectroscopy of triplet states of Rb 2 by femtosecond pump-probe 
photoionization of doped helium nanodroplets 
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The dynamics of vibrational wave packets in triplet states of rubidium dimers (Rb2) formed on 
helium nanodroplets are studied using femtosecond pump-probe photoionization spectroscopy. Due 
to fast desorption of the excited Rb2 molecules off the droplets and due to their low internal tem- 
perature, wave packet oscillations can be followed up to very long pump-probe delay times > 1.5 ns. 
In the first excited triplet state (1) 3 E^~, full and fractional revivals are observed with high contrast. 
Fourier analysis provides high-resolution vibrational spectra which are in excellent agreement with 
ab initio calculations. 
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I. INTRODUCTION 

The dynamics of vibrational wave packets (WP) has 
been investigated for a large variety of molecular systems 
ranging from simple diatomics to biomolecules [3, 0, Q • 
Alkali diatomic molecules were among the first molecules 
to be studied due to their strong electronic transitions in 
the visible range of the spectrum. Low ionization po- 
tentials make alkali dimers accessible to photoionization 
(PI) using comparatively low laser intensities. Besides, 
potential energy curves can be calculated with high preci- 
sion, thus facilitating the interpretation of spectroscopic 
data (see e. g. Q and references therein). 

A number of interesting phenomena have been studied 
on these simple molecules, e. g. WP propagation in spin- 
orbit-coupled states I5L \&, the competition of different 
ionization pathways [Tj, |8j, fractional revivals of vibra- 
tional WPs [t| [l(| , and isotope-selective ionization [ll[ . 
Detailed insights into the vibrational dynamics have been 
obtained by applying new experimental techniques such 
as photo-electron spectroscopy fl2l|. and optimal control 
schemes using shaped laser pulses [H, [l4| • 

Besides, valuable spectroscopic information can be 
gained from time-resolved experiments [H, @, Il5| . The 
extraction of high resolution spectra and potential en- 
ergy curves from ultrafast laser experiments has been dis- 
cussed in detail by Gruebele and Zewail, using molecular 
iodine I2 as a model system (l6| . 

Alkali dimers have attracted new interest in the past 
years due to the recent advances in the formation 
of ultracold molecules out of ultracold atomic ensem- 
bles by means of Feshbach resonances [l?} and photo- 
association [TH, [l!| . These studies require the knowledge 
of molecular spectra with great precision [2(| HH • How- 
ever, using conventional spectroscopy of molecules in the 
gas-phase which are formed in the singlet ground state, 
triplet transitions are more difficult to interrogate. 

Pump-probe spectroscopy of alkali dimers attached 
to He nanodroplets has been reported before by our 
group [13, HH]. Helium (He) nanodroplets are widely 
used as a nearly ideal cryogenic matrix for spectroscopy 



of embedded atoms, molecules, and clusters |24j,[25|. Al- 
kali atoms and molecules represent a particular class of 
dopant particles characterized by their extremely weak 
binding to He droplets with binding energies of about 
10 K [261 ] . From both theory and experiment it is known 
that alkali dimers reside in bubble-like structures on the 
surface of He droplets [H, [2?], [H, [29[ . Upon laser exci- 
tation, alkali atoms and molecules mostly desorb off the 
droplet surface as a consequence of repulsive interactions 
with the He surroundings induce d by the change of elec- 
tronic and vibrational states [3(1 [U [HJ . 

The influence of the He droplet on the vibrational 
dynamics of attached molecules has been studied using 
potassium dimers (K2) in singlet states 22]. Transient 
modulations of both amplitudes and frequencies of vibra- 
tional frequency components were observed, from which 
the time constant for the desorption dynamics was esti- 
mated to range between 3 and 8 ps. WP propagation in 
triplet states was observed for the first time using sodium 
dimers (Na2) attached to He nanodroplets (23|. However, 
the triplet-state dynamics did not reveal any influence of 
the He matrix on the WP dynamics. This may be due to 
considerably shorter desorption times of alkali dimers in 
triplet states as a result of more efficient coupling to the 
surrounding He as compared to singlet states. Bovino et 
al. have calculated enhanced vibrational quenching cross 
sections by orders of magnitude for vibrational states of 
the lowest triplet state a 3 £+ of Li2 as compared to the 
singlet ground state [33]. Furthermore, recent studies of 
the cw excitation spectra of Rb2 on He droplets provides 
indications that triplet dimers are oriented parallel to 
the He surface and thus may couple more efficiently to 
surface excitations of the He droplet [3^ ]. 

This work presents a systematic study of the pump- 
probe dynamics of Rb2 formed on He nanodroplets in 
triplet states in the entire tuning range of the Tita- 
nium:Sapphire (Ti:Sa) laser between A = 700 nm and 
A = 1030 nm. Only weak perturbations of the vibra- 
tional dynamics by the He matrix are observed at short 
delay times t < 10 ps associated with the desorption dy- 
namics. Thus, the dynamics probed in our experiment 
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FIG. 1: Selected triplet potential energy curves of neutral 
Rb2 and of the Rb J ionic ground state relevant to the present 
study. The arrows in (a) and (b) indicate two different path- 
ways leading to ionization: (a) Creation of a wave packet in 
the first excited state followed by resonant two-photon ioniza- 
tion; (b) Creation of a wave packet in the lowest triplet state 
by impulsive Raman scattering followed by three-photon ion- 
ization. 



corresponds mostly to the one of free Rb2 molecules in 
the gas-pase. Due to the low internal temperature of 
Rb2 formed on He droplets, exceptionally long coherence 
times of the vibrational WP oscillations are observed. 
This allows to obtain high-resolution Fourier spectra of 
the vibrational frequencies of triplet states of Rb2 , which 
may be of considerable interest in the context of the for- 
mation of ultracold Rb2 molecules 0, |35[ . 

The paper is organized as follows: Sec. HT1 gives a short 
introduction into the theory of vibrational WPs in an- 
harmonic potentials. Sec. IIIII outlines the experimen- 
tal arrangement. The discussion of experimental results 
begins with the presentation of PI and beam depletion 
spectra in Sec. IIV1 It is followed by a more detailed dis- 
cussion of the temporal and spectral data measured in 
different laser wavelength ranges corresponding to three 
different electronic transitions (Sec. lVirVIIIj) . A quantita- 
tive analysis of the pump-probe spectra measured around 
A = 1000 nm including the comparison with theoretical 
potential curves is given in Sec. IVII 



II. VIBRATIONAL WAVE PACKET DYNAMICS 

The dynamics of vibrational WPs is studied using the 
pump-probe technique. The principle of operation is the 



following. The first pump pulse creates a vibrational WP, 
tp{x,t) — '^2 n a n \n)&cp—iE n t/h, i.e. a coherent super- 
position of vibrational states of a bound electronic state, 
as illustrated in Fig. [T] for the case of Rb 2 . Here, \n) de- 
notes the set of vibrational eigenstates with energies E n , 
and a n are constant coefficients. A second fs laser pulse 
probes the WP through excitation to a higher electronic 
state |/) which is detected with time-independent meth- 
ods, i. e. by measuring the yield of photo-ions. In the 
scheme shown in Fig. [T] (a) the probe step consists of a 
resonant two-photon transition towards the ionic ground 
state 2 £+ via the 4 3 S+-state of Rb 2 . The ion signal S(t) 
is detected as a function of the time delay t between the 
two pulses. 

According to the Franck-Condon principle, the WP is 
detected near one or a few internuclear separations (Con- 
don points), at which the energy difference between up- 
per and lower potential curves match the laser energy. As 
the WP swings back and forth inside the potential well, 
excitation to the ionic state becomes maximum each time 
the WP passes a Condon point. This causes the detec- 
tor signal to be modulated periodically. In case there is 
only one Condon point which coincides with one classi- 
cal turning point, the transition probability to the ionic 
state is particularly high due to the favorable overlap be- 
tween wave functions in lower and upper states (Franck- 
Condon factor, FCF). This condition is fulfilled in the 
scheme shown in Fig. [T] (a). In this case, the period 
of the ion signal approximately equals the classical fre- 
quency of the molecular vibration, oj c . The measured sig- 
nal S(t) = \P{t)\ 2 is given by P{t) = {f\ft ■ %(z, t)) = 
b n exp — iE n t/h, where (/ represents the final ionic 
state, p, the transition dipole moment, E the electric field 
vector of the probe laser, and b n = a n {f\fl- E\n). There- 
fore the signal S(t) = J2 n , m b n b mCOs [(E n - E m )t/h] is 
composed of beat frequencies between all pairs of energy 
levels that make up the WP. The lowest frequencies are 
given by the beats between neighboring vibrational levels 
\n) and \m), where n — mil. These frequencies approx- 
imately match cu c . Beat notes between further-separated 
levels, n = m±£, where I > 1, generate higher harmonics 

ofw c [nj. 

WP propagation in Morse-like anharmonic poten- 
tials has been studied theoretically by Vetchinkin et 
al. [13, HE H3 ■ The energy levels of a Morse potential 
V(R) = E D (1 - exp [-a(R - R e ))f are given by E v = 
Hoj[(v + l/2)-x e (v + l/2) 2 ]. Here, u = 4E D x e /h = 

ci\J2Ed / ji stands for the harmonic oscillator frequency 
where Ed is the dissociation energy, x denotes the anhar- 
monicity constant, a is the range-parameter of the Morse 
potential, and /i = mRb/2 is the reduced mass of the Rb2 
molecule. 

In anharmonic potentials the classical-like behavior de- 
cays with time due to dispersion of the initially created 
WP. In the Morse potential, however, the classical be- 
havior partly and fully recurs at certain times, known 
as full and fractional revivals, respectively. The time 
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FIG. 2: Schematic overview of the experimental setup. The 
upper part shows the optical system to produce femtosecond 
pulse pairs, the lower part shows the He droplet beam line. 



at which the WP is fully reconstituted with its original 
phase is given by T rev = 2h/(AE n — AE m ) — 2Tt/(ojx e ) 
for any m = n + 1, where AE„ — E n — E n -\. Frac- 
tional revivals occur at times t/T rev = p/q, where 
p/q is an irreducible fraction of integers [38|. It can 
be shown that the full revivals of the WP that oc- 
cur at t = T rev /2, T rev , 3T rev /2, . . . result from lowest- 
frequency beats [i — 1), while half-period revivals at 
t = T rev /A, T rev /2, 3T re „/4, . . . result from second har- 
monic beats {1 = 2), etc. [10| . 



III. EXPERIMENTAL 

The experimental setup is composed of two main parts, 
as schematically represented in Fig. [21 A molecular beam 
line produces He droplets doped with Rb 2 and an optical 
system supplies pairs of fs pulses with variable time de- 
lay. The beam line consists of four differentially pumped 
vacuum chambers. In the source chamber a continuous 
beam of He droplets is formed in a supersonic expansion 
of high purity 4 He gas out of a cold nozzle (T w 17 K, 
diameter d = 5 fim) at high pressure (p w 50 bar) . Un- 
der these conditions He droplets contain around 8000 He 
atoms [24[. The source chamber is pumped by three 
turbo pumps (1000 1/s each) backed by a roots and rotary 
vane pump. 

The beam enters the oven chamber through a commer- 
cial skimmer with a diameter of 0.4 mm to pick up two 
Rb atoms out of atomic Rb vapor inside a pickup cell. 
The number of collisions of the droplets with free dopant 
atoms inside the vapor cell is controlled by the oven tem- 
perature. Given the flight distance of the droplets of 
1 cm inside the vapor cell the Rb reservoir is heated to a 
temperature T «380K to achieve highest probability for 



pick-up of two dopant atoms per droplet. Due to their 
high mobility on the superfluid He droplets the pickcd- 
up atoms bind together to form dimers which are weakly 
bound in dimples on the droplet surface. Collisional as 
well as binding energy is dissipated by the He droplet 
through evaporation of He atoms, which may also cause 
the desorption of the alkali dimers from the droplets. 
Since the amount of internal energy released upon for- 
mation of ground state (X 1 S+) dimers greatly exceeds 
the one released upon formation of dimers in the low- 
est triplet state a 3 S+, the later has a higher chance to 
remain attached to the droplet. This leads to an en- 
richment of the droplet beam with high-spin dimers and 
clusters compared to covalently bound systems [H, • 
The alkali-He droplet complex eventually equilibrates at 
the terminal temperature of pure He droplets of 380 mK. 
Thus, only the lowest vibrational state, v = 0, and a few 
rotational states (J < 10) are populated which provides 
well-defined starting conditions for pump-probe experi- 
ments. 

In the third vacuum chamber the doped droplet beam 
is intersected by the laser beam at right angle inside 
the detection volume of a commercial quadrupole mass 
spectrometer. The fs laser pulses are generated by a 
commercial mode-locked Ti:Sa laser (Chameleon, Coher- 
ent) at 80 MHz repetition rate and at an average power 
output between 1.5 and 3.2 W depending on the wave- 
length. The pulses have a duration of «160 fs and a 
full spectral bandwidth at half maximum (FWHM) of 
Aujf a Y HM ~ 80 cm" 1 . The beam is sent through a 
Mach-Zehnder interferometer to split one pulse into two 
time-delayed pulses with equal intensity and phase. The 
time delay between the pulses is controlled by a com- 
mercial translation stage. In the experiments time incre- 
ments of 30 fs are used. Behind the interferometer the 
laser beams are recombined with parallel polarizations 
and focussed into the doped He droplet beam using a 
lens of 150 mm focal length. The l/e 2 -beam diameter in 
the focus is 41 /jm which leads to a pulse peak intensity 
between 0.5 and 4 GW/cm 2 depending on the wavelength 
of the laser. The fourth chamber contains a Langmuir- 
Taylor surface ionization detector that is used to check 
the do ping level and to perform beam depletion measure- 
ments djim. 



IV. PHOTOIONIZATION SPECTRA 

In a first series of measurements, PI spectra are 
recorded with mass selection of the photo-ions using only 
one fs laser pulse at varying laser wavelength. PI spec- 
tra of Rb atoms, Rb 2 dimers and Rb3 trimers attached 
to He nanodroplets are shown as solid lines in Fig. [3] 
(a)-(c). For comparison, we have recorded PI spectra 
of Rb and Rb 2 in the gas-phase (dashed lines in Fig. [3] 
(a) and (b)). Note that all spectra shown in Fig. [3] 
are raw data that do not account for varying laser in- 
tensity (Fig. [3] (d)). The gas-phase spectrum of Rb 
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FIG. 3: Comparison of fs photo-ionization spectra of Rb 
monomers, dimers and trimers attached to helium nan- 
odroplets ((a)-(c)) with beam-depletion spectroscopy (d). 
The dashed lines in (a) and (b) are spectra of Rb and Rb2, 
respectively, in the gas-phase. The shaded area in (c) indi- 
cates the absorption band (2) 4 E' <— (l) 4 ^ °f Rb3 measured 
by cw-spectroscopy [42l . |43 |. The spectra are not corrected 
for varying laser intensity (e). 



is dominated by a broad maximum around the D2-line 
(5P 3/2 <-5Si/ 2 ) at 780 nm. The fs PI signal at that 
wavelength is additionally enhanced by the transition 
5D3/2[5/2] 4 ~5P 3 / 2 which is blue-detuned from the D2- 
line by only 67.4 cm' 1 [70.4 cm -1 ] and therefore simulta- 
neously near-resonant given the bandwidth of the fs laser. 
A small peak around the transition 5D 3 / 2 *—5Pi/2 is a l so 
visible, reflecting resonance-enhanced three-photon ion- 
ization (RE3PI) via the 5Pi/2 and 5D 3 / 2 states. Interest- 
ingly, the Dl-transition (5Pi/2 «— 5S1/2) is n °t observed 
as increased PI signal, presumably due to missing further 
resonance conditions for the transition up into the ionic 
state. 



For PI of Rb attached to the surface of He nan- 
odroplets, however, the situation is different. Here, a 
broad feature red-shifted to the D 1-line is measured and 
the peak close to the 5D 3 / 2 ^5Pi/2-transition is also en- 
hanced. The difference between gas-phase and He droplet 
spectra is attributed to the interaction between the He 
environment and the Rb dopant which leads to shift- 
ing and broadening of atomic transitions (3^|. These 
matrix-induced perturbations are expected to be partic- 
ularly pronounced for highly excited electronic states of 
Rb due to their extended orbitals and higher polarizabil- 
ities. Thus, apparently resonant-multi-photon ionization 
conditions are met at the wavelength of the Dl-line of 
Rb attached to He droplets. 

The PI spectra recorded at the mass of Rb2 are de- 
picted in Fig. [3] (b). In order to assign the observed 
features of the spectrum of Rb2 on He droplets (red line) 
we have also recorded the PI spectrum of gas-phase Rb2 
dimers formed inside the doping-cell heated to 380 K 
by three-body collisions (dashed line). Gas-phase Rb2 
molecules are expected to occupy ro-vibrational states 
of the electronic ground state X : E+ according to their 
internal temperature. In the entire tuning range of the 
laser we observe one broad feature which is assigned to 
the A 1 ^ «— X 1 S+-system and is analyzed in Sec. IVIII 
In contrast, the He droplet spectrum reveals 3 broad 
peaks. This qualitative discrepancy between gas-phase 
and He droplet spectra reflects the different formation 
mechanisms of Rb2 molecules out of single Rb atoms. In 
contrast to strongly bound Rb2 formed in the hot vapor 
of the doping cell, the association of alkali dimers on the 
surface of He droplets favors weakly-bound states, i. e. 
the lowest triplet state a 3 E+. We therefore assign the 
additional features showing up in the droplet spectrum 
around 1000 nm and around 750 nm to the triplet sys- 
tems (1) 3 S+ <- a 3 S+ and (l) 3 n g <- a 3 E+, as detailed 
m sections fVl and IVlTTl 

In order to assign the center feature of the Rb2 spec- 
trum around 850 nm we have also measured the PI spec- 
trum of Rb 3 trimers formed on He droplets, shown in 
Fig. ||c). Fragmentation of Rb^ -^Rb^+Rb upon PI 
may cause a crosstalk of the Rb^ ion signal onto the 
RbJ signal and thus perturb the Rb 2 PI spectrum. Rb 3 
photo-ions are detected nearly in the entire scan range of 
the laser. High-spin Rb 3 trimers formed on He droplets 
have been studied using cw spectroscopy [H, |43[. In 
the spectral range of interest, one electronic transition 
was identified ((2) 4 _E' <— (1) 4 A 2 ) in the range between 
11500 and 11700 cm -1 , indicated by the shaded area in 
Fig. [3] (c). Indeed, the He droplet PI spectra measured 
on Rb2 and Rb 3 masses in the range between 820 and 
900 nm show qualitative resemblance. This suggests that 
fragmentation of Rb 3 contributes to the Rb2 spectrum. 
Besides, in this spectral range the Rb 2 -He droplet spec- 
trum partly coincides with the system A <— X measured 
in the gas-phase. Therefore we conclude that presum- 
ably both the singlet transition A <— X of Rb2 as well 
as fragmentation of Rb 3 contribute to the PI spectrum 
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measured at the mass of Rb2. 

Finally, Fig.[3](d) shows the beam depletion (BD) spec- 
trum recorded using Langmuir- Taylor surface ionization 
detection. This measurement is performed using the un- 
focused fs laser beam to maximize the spatial overlap of 
the laser with the droplet beam. Thus, laser intensity is 
lower by about a factor 10 -3 as compared to the PI mea- 
surements mentioned above. We take as depletion signal 
the relative drop of surface ionization count rate when 
shining in the fs laser beam in proportion to the signal 
measured with blocked laser beam. Nearly all of the spec- 
tral features of the PI spectra recorded on masses of Rb, 
Rb2 and Rb3 are qualitatively reproduced. This clearly 
shows that the PI spectra stem from rubidium oligomers 
attached to He nanodroplets. The fact that the exact 
shapes and relative intensities of individual peaks differ 
from the ones seen in the PI spectra is mainly due to 
the much lower laser intensity. As a consequence, only 
the lowest lying transitions are expected to contribute to 
the signal, whereas the PI signal involves multi-photon 
resonances. Even though the detailed comparison of PI 
and BD spectra may provide new insight into He matrix- 
induced shifts and broadenings of higher lying electronic 
levels, it lies beyond the scope of this paper. In the follow- 
ing, the vibrational WP dynamics in 3 electronic states 
of Rb2 are investigated using the one-color pump-probe 
technique. Wave packet oscillations of Rb3 are also ob- 
served at laser wavelength 850 nm and will be presented 
in a forthcoming publication. 



V. EXCITATION AROUND 1000 nm 

The most pronounced WP oscillations are observed at 
a laser wavelength around A =1000 nm. Fig. Q] (a) dis- 
plays the pump-probe transient measured at A = 970 nm. 
The inset shows a detailed view of the same data in the 
delay time range between t = 190 and 250 ps. The data 
nicely resolve WP oscillations with a period T a rs 2.5 ps 
corresponding to the WP dynamics in the lowest triplet 
state a 3 S^- The oscillation amplitude A a slowly de- 
creases but shows no periodic modulation. In contrast, 
the WP oscillations recorded at A = 1006 nm, corre- 
sponding to dynamics in the first excited state (1) 3 S+, 
are strongly amplitude-modulated (Fig. Q] (b)). The pe- 
riod is T^ g k, 0.95 ps. After an initial decrease of the os- 
cillation amplitude As g in the time range between and 
50 ps the signal builds up again to form a local maximum 
at about 80 ps which corresponds to the first half-period 
fractional revival. The first full revival occurs at about 
160 ps. Up to about 10 full revivals are observed at delay 
times up to > 1.5 ns. The data shown in the inset illus- 
trate the transition from the full revival to one-third and 
one-half period fractional revivals of the WP dynamics 
in the (l) 3 E+-state. The initial phase 4>o, I006nm of the 
WP oscillation at A = 1006 nm can be obtained from the 
extrapolation to t = of the linear fit of the positions of 
signal maxima as a function of delay time between t = 
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FIG. 4: Experimental yield of Rbj photo-ions as a function of 
delay time between pump and probe laser pulses recorded at 
laser wavelength 970 nm (a) and 1006 nm (b). The insets show 
detailed views in the time range of the transition from full to 
half-period fractional revivals of the wave packet dynamics in 
the (l) 3 E+-state. 



and t = 20ps. The resulting value 4>o, ioo6nm = 0.45-27ris 
consistent with the excitation scheme depicted in Fig. [T] 
(a), which illustrates the propagation of a WP created 
at the right turning point of the (l) 3 £+-state potential 
and probed at the left turning point after 1/2 oscilla- 
tion period. In addition to the periodic modulation due 
to revivals, the oscillation amplitude A^ g is damped fol- 
lowing roughly an exponential decay behavior with de- 
cay time constants varying between 100 and 700 ps de- 
pending on the excitation laser wavelength. This slow 
damping is attributed to contributions of different ro- 
tational states to the measured transients, which have 
slightly shifted vibrational frequencies due to rotation- 
vibration coupling [16]. Since the rotational tempera- 
ture of the molecules formed on He nanodroplets is low 
compared to gas-cell or standard molecular beam exper- 
iments, though, this vibrational dephasing is very slow. 

As a consequence of the high-contrast WP oscillations 
that persist up to > 1.5 ns we may expect to obtain 
Fourier transform spectra with high resolution. Fig. [5] 
gives an overview of the Fourier spectra of the pump- 
probe transients recorded at selected laser wavelengths. 
One can clearly distinguish 2 different sets of frequen- 
cies: The spectrum at A = 970 nm is dominated by the 
fundamental frequency component uj a w 13.3 cm -1 and 
higher harmonics up to uj a j w £u> a , where £ = 2,3,4. 
As detailed in Sec. [Vll this frequency is assigned to 
WP dynamics excited by resonant impulsive stimulated 
Raman scattering (RISRS) in the lowest triplet state 
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FIG. 5: Fourier transforms of wave packet oscillations at dif- 
ferent laser wavelengths. At laser wave length A = 970 nm 
the spectrum is dominated by vibrations in the lowest triplet 
state a 3 £j , at A = 1006 nm only the excited (l) 3 E^"-state 
dynamics is present, at A = 1025 nm both signals contribute. 
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TABLE I: Experimental values for the lowest vibrational fre- 
quencies uv °f the a-state of Rb 2 with respect to u) v ii =0 . 
The corresponding rotational constants B v n are calculated 
from the fitted Morse potential. Units are cm -1 . 



a 3 Sj according to the excitation pathway sketched in 
Fig. [U (b) [7] . The fundamental frequency corresponds 
to energy level spacings between neighboring vibrational 
levels, i.e. to a = (E v /i — E v ii _i) / (he) . Higher harmon- 
ics u a ,e — (E v » — E V "-g)/(hc), where £ = 2, 3, 4, are 
present in the spectrum as a consequence of overtone 
transitions v" = I <— v" = which are allowed due to 
the anharmonic shape of the potential V a . However, their 
amplitudes are strongly suppressed compared to the one 

of LO a . 

Fig. [5] (b) shows a closeup view of the spectrum in the 
range of the fundamental frequency u> a — (£7i — Eg) /(he) 
recorded at laser wavelength 1025 nm. The high spectral 
resolution allows to separate the fundamental frequen- 
cies co a Ri 13.2 and 13.3 cm -1 of the two isotopologues 
85 Rb2 and 87 Rb2. The small peak in between at about 
13.25 cm" 1 comes from the mixed isotopologue 85 Rb 87 Rb 



which adds to the 85 Rb2 signal as a consequence of in- 
sufficient mass resolution of the QMS detector. In addi- 
tion, small peaks shifted to lower frequencies by about 
0.4 cm -1 are visible in the spectrum as beat frequencies 
between excited vibrational states, (E^— Ex)/ (he). Their 
low peak intensity directly reflects the weak thermal pop- 
ulation of excited vibrational states v" > due to the low 
temperature of the molecules formed on He nanodroplets. 
The peak positions u> v " for 1 < v" < 5 with respect to 
luq are determined by fitting Lorentzian model curves to 
the data and are summarized in Table U The uncertainty 
associated with each value of ov is estimated from the 
precision of the fit results to be < 0.003 cm -1 . In addi- 
tion to the level frequencies, Table Q] lists the rotational 
constants calculated using the Morse potential extracted 
from the vibrational spectrum (cf. Sec. IVIj) . 

It is important to note the high spectral resolution ob- 
tained using these fs pump-probe measurements. E. g., 
the peaks in Fig. [5] (b) have full widths at half max- 
imum (FWHM) Alj^ whm w 0.03cm -1 (» 900 MHz). 
In particular, the resolution is higher or comparable 
with rotational, fine structure, and hyperfme splittings 
in the a-state. However, besides the vibrational spac- 
ings no additional beat frequencies are observed. The 
fact that no rovibrational cross-beatings of the type 



v',J'=J"+l- 



-E v /-i 



j'=j' 



■-i)/(hc) are observed 



may have two reasons. Since beatings between rovi- 
brational states with equal J', u> V ib = (-Ed',j'=j"±i — 
E v r-i_jr = j>i±i)/(hc) are statistically favored with respect 
to the cross-beatings w r cwi6> summing over all frequency 
contributions (assuming thermally populated J'-states) 
effectively reproduces the pure vibrational spectrum on 
top of a broad background. Additionally, the initial 
coupling to the He droplet prior to desorption may in- 
duce rotational decoherence and thereby quench the rovi- 
brational beating signals. The fact that the pure rota- 
tional dynamics with frequencies uj ro t — (E v i ,j'=j"+i — 
E v < ,j>=j» -i) / (he) — 2B V (J' + 1) is completely missing 
in all our measured transients for delay times > 10 ps 
supports the latter interpretation. 

In both cases, the contribution of different rotational 
states to the signal only becomes observable indirectly 
through the rotation-vibration coupling, which leads 
to slightly different vibrational frequencies for differ- 
ent J'-states. Summing over the J'-state contributions 
then leads to a slight broadening and shifting of the 
band origins by approximately Aw « —a e k^T/(B v hc), 
where a e denotes the rotation-vibration coupling con- 
stant and T stands for the rotational temperature [l6l ]. 
From ab initio potentials (Sec. IVI[) we find the value 
a a e w 1.6 • 10 cm" 1 , which would lead to a system- 
atic shift of the determined vibrational frequencies by 
Aoj a « —0.004 cm -1 when assuming full thermalization 
to the He droplet temperature of T = 0.4 K. 

At laser wavelength A = 1006 nm, the Fourier spec- 
trum in Fig. [5] (center panel) shows the fundamental fre- 
quency ui-zg = (E v > — E v i-i)/(hc) ~ 35 cm -1 and higher 
harmonics we 9 ^ = (E v > — E v >—e)/ (he) up to £ = 3. These 
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FIG. 6: Detailed views of the Fourier spectra around the fun- 
damental frequencies of wave packet dynamics in the first ex- 
cited triplet state (1) 3 E^" (a) and in the lowest triplet state 
a 3 E^ (b). In panels (a), different colors represent measure- 
ments at different laser wave lengths. The vibrational fre- 
quencies of the isotopologues 85 Rb 2 and 87 Rb 2 are consider- 
ably shifted. 



frequencies reflect the wave packet propagation in the 
first excited triplet state of Rb2, (1) 3 E+, excited accord- 
ing to the scheme shown in Fig. Q] (a). At A = 1025 nm, 
both sets of frequencies are present with roughly equal 
intensities. 

The detailed view of the Fourier spectrum around us 9 , 
shown in Fig. [5] (a), reveals progressions of vibrational 
frequency differences for the two isotopologues 85 Rb2 
and 87 Rb2- Different peak colors indicate measurements 
using selected laser wavelengths leading to the excita- 
tion of different sets of vibrational states v' according to 
energy conservation and Franck-Condon principle. All 
individual spectra are normalized to the height of the 
dominant peak. From these Fourier spectra frequency 
differences u> v > = (E v > — E v i_\)/(Jic) between neighbor- 
ing levels v 1 — (0, 1) up to v' — (12, 13) can be de- 
termined with high precision. Table |TT] summarizes the 
values of u) v i obtained from Lorcntzian fits with an un- 
certainty < 0.005 cm -1 . In addition, rotational con- 
stants calculated from Morse potentials obtained from 
the measured data are listed (Sec. IVI[) . The systematic 
shift due to rotation-vibration coupling is estimated to 
Awsg ~ — 0.0015 cm -1 , given the rotation- vibration cou- 
pling constant of 9 « 5.8 ■ 10~ 5 cm" 1 and the rotational 
temperature T = 0.4 K. 

The temporal behavior of individual spectral compo- 





8b Rb 2 


8V Rb 2 


/ 

V 


AiU v i 


B v r 


Aov 


B v > 







0.01486 




0.01452 


(0,1) 


36.382 


0.0148 


35.966 


0.01446 


(1,2) 


36.179 


0.01474 


35.765 


0.0144 


(2,3) 


35.972 


0.01468 


35.565 


0.01434 


(3,4) 


35.763 


0.01462 


35.361 


0.01429 


(4,5) 


35.555 


0.01456 


35.156 


0.01423 


(5,6) 


35.345 


0.0145 


34.945 


0.01417 


(6,7) 


35.162 


0.01444 


34.764 


0.01411 


(7,8) 


34.953 


0.01438 


34.558 


0.01405 


(8,9) 


34.737 


0.01432 


34.350 


0.01400 


(9,10) 


34.523 


0.01426 


34.139 


0.01394 


(10,11) 


34.332 


0.0142 


33.957 


0.01388 


(11,12) 


34.119 


0.01414 


33.745 


0.01382 


(12,13) 


33.910 


0.01408 


33.533 


0.01376 


(13,14) 


33.714 


0.01402 


33.331 


0.01370 



TABLE II: Experimental level spacings ui v r = {E v i — 
E v i_ 1 )/(hc) and calculated rotational constants B v i of the 
lowest vibrations v' of the first excited state (1) 3 E+ of Rb 2 . 



nents can be visualized by means of the spectrogram rep- 
resentation Fig. [7] shows the same data as Fig. [5j 
analyzed by Fourier transformation inside a 10 ps-time 
window of Gaussian shape which slides across the entire 
data set. The vertical dark lines in the top panel repre- 
sent constant frequencies uj a ,e as a function of delay time. 
The fading away of signal intensity at long delay times 
> 100 ps reflects damping of the WP motion, as men- 
tioned above. Clearly, amplitudes of the higher harmon- 
ics u> a ,t) where £ > 1, fall off faster as compared to the 
fundamental frequency component Lu a . At A = 1006 nm, 
full and fractional revivals of the wave packet dynamics 
nicely show up as periodic intensity maxima as a function 
of delay time. In this representation, half-period and one- 
third-period revivals occurring at t 1 / 2 = k ■ T rev /A and 
at t 



1/3 



k ■ T rev /6, where k = 1,2,3, 



respectively, 



are clearly discerned as second and third harmonic beats 
(see Sec. IT]). Again, WP coherence fades away increas- 
ingly quickly for the higher harmonics to-zg^. 

In contrast to the pronounced revival structure of the 
(l) 3 £+-state dynamics, the a-state dynamics features 
only very weak revivals at integer multiples of about 80 ps 
(Fig. [4]). This is due to the different distributions of pop- 
ulation in vibrational level manifolds v" and v'. Since 
in the a-state predominantly v" = is populated by the 
pump pulse and to a much lesser extent v" = 1 and higher 
v"-states, the dynamics primarily results from the beat- 
ing between the two levels v" = and v" — 1. This 
causes WP oscillations with only weak revival structure. 
In contrast, in the (l) 3 £+-state a number of u'-levels are 
roughly equally populated which gives rise to high con- 
trast revival features. 

It is interesting to note that even when the pump- 
probe-transient is entirely dominated by WP-dynamics 
in the a-state (cf. Fig. [7] (a)), there are still contribu- 
tions from (l) 3 S+-state dynamics that rapidly disappear 
at delay times t > 10 ps (signal at oj^ g w 33 cm -1 ). 
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FIG. 7: Spectrogram representation of the vibrational wave 
packet dynamics of Rb2 at selected laser wavelengths. Con- 
tinuous vertical bands correspond to monotonous oscillation 
amplitudes, dotted bands reflect revivals of the wave packet 
motion. 



We take this as an indication that, while WP-dynamics 
at t > 10 ps reflects unperturbed gas-phase vibrations, 
frequencies and amplitudes of WP-oscillations at short 
delay times are affected by He matrix interactions. In 
particular, the unexpectedly prominent a-state dynamics 
associated with RISRS-excitation may be due to strong 
quenching and decoherence of vibrations in the excited 
(l) 3 S+-state by the He droplet as compared to damping 
of a-state vibrations. The effect of the He environment 
on the rovibrational dynamics of different ground and ex- 
cited states of Rb2 and other molecules will be further 
studied both experimentally and theoretically. 



VI. ANALYSIS OF FOURIER SPECTRA 

In order to compare the experimental results with theo- 
retical predictions, the experimental data are fitted to the 
energy terms of the Morse potential (see Sec. IJ). Thus, 




33.0-1 , 1 , 1 , 1 , 1 , 1 , 1 , ,- 

2 4 6 8 10 12 14 

Vibration v' in (if z* -state 



FIG. 8: Linear fits of the experimental frequency differences 
between adjacent vibrational levels of the (l) 3 E+-state of 
85 Rb 2 and 87 Rb 2 . 



State 


Rb 2 


E D /cn\' 1 


ABo/cm -1 


a/A- 1 


Aa/A- 1 


a 


85 


247.90 


0.40 


0.69044 


0.0006 




87 


244.3 


3.2 


0.6963 


0.005 


(1) 3 S+ 


85 


3246.42 


4.3 


0.50954 


0.0003 




87 


3238.015 


0.94 


0.51023 


0.00007 



TABLE III: Morse potential parameters and absolute uncer- 
tainties obtained from fitting the experimental frequencies by 
the Morse term energies. 



by determining lu and Ed from the experimental data 
a Morse potential can be specified as approximation to 
the real molecular potential curve. Only the position of 
the potential well, i? e , remains undefined and will be ad- 
justed by comparing with ab initio potentials published 
as supplementary material in Ref. [201 ] . Energy differ- 
ences between neighboring levels are given by the linear 
dependence AE = E v — E v _x = hu(l — vuj/(2E d )). 
The linear fits of the experimental values of ws 9 for the 
(l) 3 S+-state of the two homonuclear isotopologues of 
Rb2 are shown in the Birge-Sponer plot (Fig. [H]). The 
error bars are barely visible on the given frequency scale. 
The agreement between the experimental data and the 
linear model is excellent in the energy range of E v r for 
< v' = < 13. The same analysis is performed on 
the data for the a-state. The resulting Morse potential 
parameters are summarized in Table Hill The associated 
potential curves for the (1) 3 £+ and a-states of 85 Rb2 
and 87 Rb2 are plotted in Fig. [9] (a) and (b) as dashed 
and dotted lines, respectively. The good agreement of 
the curves for the two isotopologues highlights the high 
level of precision of the determination of vibrational fre- 
quencies from our data. 

In addition to the Morse potentials, Fig. [TJ] shows the 
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FIG. 9: Morse potentials of the first excited triplet state 
(1) 3 E+ (a) and of the lowest triplet state a 3 Ej (b) derived 
from the measured vibrational frequencies compared to ab 
initio calculations. The insets show magnified views of the 
long-range parts of the potentials. 



ab initio potential curves as black solid lines. As rec- 
ommended in Ref. [2Cj, a short-range effective poten- 
tial has to be added to the computed potential curves. 
Therefore, the ab initio (l) 3 £+-state potential is modi- 
fied by adding an additional repulsive potential V sr (R) — 
1.613181e- 3 - 5i lR 8 - 5 according to Ref. 0. In the same 
way the a-state potential curve is corrected by adding 
a short-range potential, which is tabulated in Ref. [i^ . 
As expected, close to the dissociation limit the Morse 
potentials significantly deviate from the shape of the ab 
initio potential. However, the depths of the potential 
wells Ed of both a- and (l) 3 E+-states are in remarkably 
good agreement with the ones of the ab initio potentials, 
E D,ab initio /( hc ) = 3313.7cm- 1 and E% ab initio / '{he) = 
254.9 cm -1 , considering the limited energy range cov- 
ered by the observed vibrational levels (straight horizon- 
tal lines in Fig. [5]). 

The fact that in our experiment the revival structure of 
the (l) 3 £+-dynamics is very pronounced and long-lived 
is the consequence of the good agreement of the (1) 3 S+- 
state potential curve Vs g with the Morse potential. Devi- 
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FIG. 10: Dependence on laser intensity of the first exciting 
pump pulse of the fringe contrast of a and (l) 3 S^"-state dy- 



ations from the shape of the Morse potential would cause 
the revival structure itself to dephase and the WP would 
only be fully recovered at much longer times. This effect, 
however, is not at the origin of the observed exponential 
damping of the oscillation amplitudes. 

Using the fit parameters we can now determine the 
delay time when the first full revival is expected to 
occur, t = T rev /2 = ir/(u!X e ) — AwE D / (Hu 2 ). For 
85 Rb 2 [ 87 Rb 2 ] in the (l) 3 S+-state we find the values 
T rev /2 =161. 8ps [165. lps], which are in excellent agree- 
ment with the observation. For the a-state dynamics we 
calculate revival times T rev /2 =90.6ps [90. Ops], which 
approximately match the observed values. 

Assigning the two sets of frequencies to WP dynamics 
in the lowest and first excited states is further confirmed 
by selectively measuring the amplitude of WP dynam- 
ics in a and (l) 3 E+-states as a function of varying laser 



intensity of the first exciting pump pulse, J p 



To se- 



lectively attenuate the pump pulse a variable attenua- 
tor is introduced in one of the two arms of the Mach- 
Zchnder interferometer. Great care is taken to realign 
pump and probe laser beams coaxially for every setting 
of the attenuator. The a-state dynamics is selectively 
probed by measuring the maximum pump-probe signal 
at laser wavelength A = 950 nm, at which the dynamics 
is nearly entirely determined by WP propagation in the 
a-state. The delay time t = 80 ps is chosen such that 
transient effects by the interaction of Rb 2 with the He 
environment during desorption can be safely neglected. 
The (l) 3 E+-state dynamics is probed at A = 1006 nm at 
the maximum of the first full revival (t = 160 ps). In 
order to isolate the pure effect of the pump pulse in the 
pump-probe scheme the background count rates result- 
ing from direct 3PI when applying only the pump pulse 
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FIG. 11: Comparison of the calculated Franck- Condon factors 
(FCFs) for the transition (l) 3 E+(«') <- a(v" = 0) with the 
measured photoionization (PI) spectrum, with the beam de- 
pletion (BD) spectrum and with the amplitudes of the Fourier 
transform (FT) frequency components of a and (l) 3 E^"-state 
wave packet dynamics. 



or only the probe pule are subtracted from the data. 
The resulting state selective signals are shown as sym- 
bols in Fig. [10] The solid lines represent fit curves to 
the experimental data to guide the eye. The w a -data 
are modeled by a second order polynomial to account 
for the second order process of the RISRS scheme. The 
WEg-data are modeled by a simple saturation model to 
account for saturation of the strong resonant transition 
(1) 3 E+ <— a. Clearly, the data follow qualitatively dif- 
ferent behavior which confirms the excitation schemes 
shown in Fig. [T] and thus supports the assignment of the 
observed frequency components presented in Fig. [5] to a- 
and (l) 3 £+-state vibrations. 

The correct identification of the (l) 3 £+-state WP dy- 
namics is further confirmed by comparing the one-pulse 
Pi-spectrum shown in Fig. [3] in the wavelength range 
around 1000 nm with calculated Franck-Condon factors 
(FCFs) for the transition (l) 3 £+(t/) <- a(v" = 0). 
The FCFs are computed using LeRoy's LEVEL 8.0 pro- 
gram |4j| with the ab initio potential as input. The re- 
sulting values are plotted in Fig. [11] as black diamonds, 
the experimental spectrum is shown as red solid line. The 
agreement between measurement and calculation regard- 
ing peak positions and widths is quite good, in particu- 
lar when considering additional broadening of the fs PI 
spectrum due to the spectral bandwidth of the fs-laser 
as well as saturation effects. However, it has to be as- 
sumed that this agreement relies on the accidental can- 
celation of He-matrix induced frequency shifts and the 
effect of near-resonant higher lying electronic states of 
Rb2, since the peak measured by beam depletion (BD) at 
lower laser intensity (blue dotted line) is blue shifted by 
about 100 cm -1 with respect to the PI maximum. Thus, 
the discrepancy between the line positions in the PI and 



BD spectra is attributed to enhanced ionization prob- 
ability due to higher lying resonances at certain laser 
wavelengths. 

The additional red squares and blue circles represent 
amplitudes of the Fourier components, A a and A^ g , of 
the fundamental frequencies ui a and respectively. 
The data have been corrected for varying laser intensity 
when scanning the laser wavelength. Interestingly, the 
A-£ g amplitude spectrum is shifted to smaller frequen- 
cies by about 200 cm -1 with respect to the one-pulse 
PI spectrum. This is an effect of the additional fre- 
quency selectivity imposed by the transition from the 
(l) 3 S+-state into the ionic continuum of Rb2. In fact, 
when incorporating into the calculation the FCFs for 
the transition from the (l) 3 £+-state to the intermedi- 
ate near-resonant state 4 3 E+, good agreement with the 
Asg-amplitude spectrum (blue circles) is achieved. 

The A a amplitude spectrum features two maxima in 
the wings of the one-pulse PI spectrum around A = 
1020 nm and A = 970 nm, and a minimum at the peak 
position of the PI spectrum around A = 990 nm. Again, 
this is attributed to the influence of various near-resonant 
higher-lying states on the transition to the ionic state 
driven by the probe pulse. This interpretation is backed 
by the analysis of the initial phases of the a-state oscilla- 
tions obtained from band-pass filtered pump-probe tran- 
sients. At A = 1020 nm, the initial phase </>o, io20nm — 
0.26 • 2ir significantly differs from the one obtained at 
A = 970 nm, </>o,97unm = 0.68 • 2-7T, which highlights the 
presence of two distinct detection windows. 



VII. EXCITATION AROUND 850 nm 

The signal of RbJ photo-ions measured on He 
nanodroplets features an additional maximum around 
850 nm, which roughly coincides with a broad feature 
measured in the effusive beam (Fig. [3] (b)). Since only 
strongly bound singlet molecules are expected to be 
formed in a hot effusive beam, these peaks are assigned to 
singlet transitions. This assignment is further supported 
by the comparison with calculated FCFs, as shown in 
Fig. [T2j The red dashed line represents the Frank- 
Condon factors for the transitions ^4 1 E,t AT 1 E+(u" = 
0) from the lowest vibrational state v" — of the singlet 
ground state X into vibrational states of the first excited 
singlet state A. This distribution is in reasonable agree- 
ment with the left peak of the droplet spectrum of Rb 2 . 
The slight deviation of peak positions may be due to per- 
turbations by the He matrix or due to the influence of 
higher lying resonances in the 3PI process. In addition, 
spin-orbit coupling between A and 6 3 n„-states are known 
to induce strong perturbations [47J. The black dashed 
line is obtained by summing over Franck-Condon spectra 
assuming thermally populated v" -levels at the tempera- 
ture of the Rb vapor cell (380 K). The resulting profile 
matches the measured PI spectrum of effusive Rb 2 quite 
well. The peak in the BD spectrum (blue dotted line) is 
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FIG. 12: (a) Potential energy curves of ground and first ex- FIG. 13: Comparison of calculated Franck-Condon factors 

cited singlet states of Rb2. (b) Photoionization spectra of Rl>2 for the transition (l) 3 O s <— a(v" = 0) with the measured 

in an effusive beam and of Rb2 attached to He nanodroplets photoionization and with beam depletion spectra, 
(solid lines) in comparison with the beam depletion spectrum 
(blue dotted line) and with calculated Franck-Condon factors 
(FCFs) for the transition A 1 ^ <- X X E+ (dashed lines). 

VIII. EXCITATION AROUND 750 nm 



narrower but centered at the same position with respect 
to the PI signal. The frequency shift of about 300 cm -1 
to the blue as compared to the calculated Franck-Condon 
profile is attributed to the effect of the He matrix. 

Besides the mentioned singlet transition, this part of 
the He droplet spectrum measured at the mass of Rb2 
may also result from fragmentation of Rbj}" — >Rb 2 +Rb 
and possibly of larger oligomers following PI, as sug- 
gested by the resemblance of Rb 2 and Rb 3 spectra (Fig. [3] 
(b) and (c)). In particular, at about 840 nm both Rb 2 
and Rb 3 droplet spectra feature coinciding local max- 
ima. Unfortunately, no WP oscillations in the pump- 
probe transients of Rb 2 are observed in this wavelength 
range. However, around 850 nm, in the region of the tran- 
sition (2) 4 E' <— (1) 4 A' 2 of Rb3 in the quartet manifold, 
we do observe WP oscillations with low contrast when 
setting the ion detector to the mass of Rb3. These data 
will be presented in a forthcoming publication. The fact 
that these oscillations are not observed on the Rb 2 -mass 
indicates that the contribution of Rb3 fragmentation to 
the Rb 2 -signal is smaller than about 20 % of the signal 
intensity. 



The Rb 2 PI spectrum displays an additional, slightly 
structured maximum around 750 nm, shown in Fig. [T3] 
(b) as red solid line. The BD spectrum (blue dotted 
line) in this spectral range significantly differs from the 
PI spectrum by featuring a much narrower peak. The 
lack of effusive signal at that wavelength already points 
at transitions in triplet molecules. Indeed, the transi- 
tion from the lowest triplet state a to the second excited 
state (l) 3 n g is found to fall into this wavelength region. 
Since the (l) 3 n s state is significantly perturbed by first- 
order spin-orbit coupling, fine-structure states 0", l g , 
and 2 g have to be considered. The corresponding po- 
tential curves are depicted in Fig. [T3] (a) including the 
vibrational levels probed in the experiment. The FCFs 
for transitions from a(v" = 0) into these states are shown 
in Fig. [13] (b) in comparison with the measured PI and 
BD spectra. While the BD spectrum agrees quite well 
with the Franck-Condon spectrum when accounting for 
the laser bandwidth, both width and fine-structure of 
the PI spectrum are not correctly reproduced. The fact 
that the line intensities of the low-frequency part of the 
spectrum (< 13400 cm" 1 ) varies much stronger as a func- 
tion of laser intensity than the high-frequency part indi- 
cates that near-resonant transitions at higher excitations 
of Rb 2 perturb the spectrum particularly strongly in the 
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FIG. 14: Fourier spectrum of the wave packet dynamics in 
spin-orbit states of (l) 3 n g . 



vicinity of the (l) 3 ILj <— a-transition. 

This transition has been studied before using cw exci- 
tation and detection of laser-induced fluorescence as well 
as magnetic circular dichroism [34| . The cw excitation 
spectrum featured 2 broad maxima around 735 nm and 
744 nm, quite similar to our BD spectrum. The inter- 
pretation by Aubock et al. was based on a molecular 
Hamiltonian containing spin-orbit coupling and a "crys- 
tal field" interaction with the surface of the He droplet. 

Pump-probe transients around 750 nm are modulated 
by WP oscillations, though the fringe contrast is much 
lower than the one observed around 1000 nm. The 
Fourier spectrum of the delay scan recorded at A = 
750 nm is shown in Fig. [14] The strongest frequency 
components are the fundamental ground state frequency 
uj a (0,l) « 13.3 cm -1 as well as w a (l,2) w 12.9cm -1 . 
Thus, the dominant excitation mechanism is RISRS. In 
addition, frequency components show up at 17.8cm -1 , 
20.4cm -1 , 22.4cm -1 , 24.1cm -1 , and at 25.6cm -1 . This 
Fourier spectrum is in good agreement with the expected 
transition frequencies and corresponding FCFs of the 
transitions 2 g , l g , ~ (v 1 ) <— a(v" — 0), depicted as 
stick spectra in Fig. [13] (b). Due to the accidental near- 



coincidence of vibrational level spacings in the three fine- 
structure states, however, it is not possible to ascertain 
which one predominantly determines the wave-packet dy- 
namics. 



IX. CONCLUSION 

In conclusion, we have systematically investigated fem- 
tosecond photo-ionization spectra of Rb2 attached to He 
nanodroplets in the wavelength range from A = 700 
to A = 1030 nm. Three transitions between low- lying 
triplet states ((1) 3 S+ <— a and (l) 3 n 9 <— a) and singlet 
states (A <— X) have been identified. One-color pump- 
probe transients around A = 1000 nm reveal pronounced 
modulations due to WP oscillations in the first excited 
triplet state (1) 3 £+ as well as in the lowest triplet state 
a 3 E+. Fractional revivals of the (l) 3 E+-state dynamics 
are visible with high contrast and long coherence times 
up to 1.5 ns due to the low internal temperature of the 
molecules formed on He nanodroplets. This allows to 
extract vibrational frequencies with high precision. The 
data are well fitted by Morse-like potential curves and 
are in good agreement with ab initio potentials. 

With this work we show that valuable spectroscopic 
information can be obtained from femtosecond spec- 
troscopy of doped He nanodroplets. In the future the 
details of the dynamic coupling of dopant molecules with 
the surrounding He quantum fluid will be further an- 
alyzed. In particular, the choice of dopant molecules 
with different coupling strengths to the He droplet, rang- 
ing from surface-bound alkali dimers to deeply solvated 
molecules, e. g. I2, in combination with the potential to 
compare superfluid He and normal fluid 3 He droplets 
offer great chances for studying the dynamics of complex 
quantum systems. 
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